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SUMMARY 


Classical  mechanisms  (e.g.,  collisions,  chemistry,  single  particle 
dynamics,  and  radiation)  always  contribute  to  the  evolution  of  nonequilibrium 
plasmas.-  From  the  theoretical,  experimental,  or  observational  viewpoints,  the 
study  of  classical  mechanisms  is  an  especially  important  first  step  which  defines 
the  upper  limit  to  the  time  scale  for  plasma  evolution  and  also  gives  the  minimum 
level  of  turbulence  required  to  make  a  plasma  anomalous.  In  fact,  the  development 
of  a  plasma  can  be  classified  as  being  anomalous  only  after  the  classical  behavior 
of  the  plasma  is  understood.  In  practice,  a  particular  plasma  can  consist  of 
different  components,  some  of  which  act  in  a  primarily  classical  way,  while  the 
rest  have  an  anomalous  behavior.  For  example,  in  the  tokamak  fusion-oriented 
device,  ion  heat  transport  is  close  to  being  classical  while  electron  heat  transport 
is  anomalous.  However,  in  the  bumpy  torus  fusion-oriented  device,  electron  trans¬ 
port  is  classical  while  ion  transport  is  anomalous. 

Classical  mechanisms  are  necessary,  basic  ingredients  in  theoretical  or 
numerical  representations  of  the  high-altitude  nuclear  environment. 

This  report  is  a  compilation  of  three  separate  papers,  each  of  which 
describes  and  analyzes  a  classical  mechanism  which  can  influence  the  ionospheric 
and  magnetospheric  environments  following  a  high-altitude  nuclear  burst.  The 
three  mechanisms  considered  are  (1)  charge  exchange,  (2)  radiative  recombination, 
and  (3)  nonadiabaticity. 

Chapter  1,  "The  Decay  of  Structure  by  Charge  Exchange,"  considers  struc¬ 
tural  decay  by  charge  exchange.  Particular  emphasis  is  given  to  those  structures, 
consisting  of  energetic  ions  with  kinetic  energies  on  the  order  of  10  or  even 
larger  than  100  keV,  which  may  develop  following  a  high-altitude  nuclear  burst  In 
the  Ionosphere.  Charge  exchange  affects  the  lifetime  of  these  structures  by 
transforming  Ions  Into  neutrals  whose  motion  Is  not  fixed  to  the  magnetic  field 
geometry.  The  result  Is  a  loss  rate  for  particles  which  does  not  depend  on  the 
magnetic  field  to  lowest  order.  This  loss  is  either  convective  or  diffusive 
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depending  on  whether  the  effective  neutral  mean  free  path  across  the  structure  is 
larger  or  smaller  than  the  structure  size  transverse  to  the  magnetic  field, 
respectively.  It  is  demonstrated  that  structures,  consisting  initially  of 
monatomic  oxygen  ions,  are  attenuated  by  at  least  two  orders  of  magnitude  in 
traveling  from  altitudes  above  400  kilometers  down  to  altitudes  below  approxi¬ 
mately  200  kilometers  where  strong  ionization  of  ambient  ionospheric  neutrals  and 
fluid  structuring  occurs.  Recent  fluid  simulations  suggest  that  density  vari¬ 
ations,  below  15  percent  on  a  larger  density  profile,  do  not  determine  the  long¬ 
term  outer  scale  length  of  the  power  spectral  density  and  so  the  attenuation  of 
structure  whose  scale  size  is  smaller  than  10  kilometers  may  not  impose  the  outer 
scale  length  at  lower  altitudes.  Indeed,  morphology,  consisting  of  energetic  ions 
with  pitch  angles  approaching  90  degrees  and  forming  near  the  geomagnetic  equator, 
imply  outer  scale  lengths  at  lower  altitudes  which  are  comparable  to  or  greater 
than  100  kilometers.  It  follows  that  structures  like  prompt  striations,  whose 
scale  size  is  comparable  to  or  smaller  than  the  ion-gyroradius,  do  not  specify  the 
long-term  outer  scale  lengths  at  lower  altitudes  by  propagation  from  starting 
points  at  higher  altitudes. 

Chapter  2,  "Analysis  of  Radiative  Recombination  Effects  on  Density  Profiles 
and  Their  Fourier  Spectrums,”  discusses  the  role  that  radiative  recombination  can 
play  in  the  modification  and  decay  of  cool,  dense,  ionospheric  morphology.  Fourier 
spectrums  of  density  profiles  are  functions  of  time  when  recombination  is  active. 
Initially,  Fourier  components  which  are  not  at  first  present  in  the  spectrum  grow 
relative  to  the  Fourier  component  which  alone  is  present  in  a  uniform  plasma. 
However,  eventually  they  and  the  other  Fourier  components  present  because  of  plasma 
nonuniformity  decay  relative  to  the  uniform  Fourier  component.  Because  the  uniform 
Fourier  component,  which  is  proportional  to  total  density,  always  decays  because  of 
recombination,  it  is  evident  that  for  long  times  the  plasma  becomes  uniform  on  a 
faster  time  scale  than  total  ionization  decay. 

Chapter  3,  "Nonadiabaticity  and  the  High-Altitude  Nuclear  Environment," 
points  out  that  nonadiabatic  changes  in  magnetic  moment  may  play  a  substantial  role 
in  the  evolution  of  the  ionosphere  and  magnetosphere  following  a  high-altitude 
nuclear  burst.  This  is  because  the  magnetic  moment  is  only  an  adiabatic  constant 
of  motion  in  magnetic  field  configurations  having  sufficiently  gentle  magnetic 
field  gradients  and  slow  temporal  variations.  However,  spatial  or  temporal  varia¬ 
tions  of  the  magnetic  field  of  sufficient  size  are  not  consistent  with  the 
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constancy  of  the  magnetic  moment.  This  nonadiabaticity  can  be  a  factor  in  the 
ionospheric  and  magnetospheric  environments,  following  a  high-altitude  nuclear 
burst,  in  two  ways.  First,  a  high-altitude  nuclear  burst  can  push  aside  the  ambi¬ 
ent  geomagnetic  field.  A  new  magnetic  field  configuration  is  created  with  sub¬ 
stantial  magnetic  field  gradients  and  local  extremums  in  magnetic  field  strength 
where  nonadiabatic  jumps  in  the  magnetic  moment  can  occur.  Because  these  jumps 
depend  on  both  pitch  angle  and  gyrophase,  nonadiabaticity  supports  the  premise  of 
outer  scale  length  being  at  least  as  large  as  the  gyroradius  of  energetic  ions. 
Second,  ions  injected  into  the  magnetosphere  following  a  high-altitude  nuclear 
burst  can  be  confined  on  a  specific  magnetic  flux  surface  in  the  magnetosphere  only 
if  their  kinetic  energy  is  below  a  bound  determined  by  nonadiabatic  considerations. 
Hence,  nonadiabatic  effects  can  influence  the  coupled  ionospheric  and  magneto- 
spheric  current  system. 
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1.  THE  DECAY  OF  STRUCTURE  BY  CHARGE  EXCHANGE 


1.1  INTRODUCTION 

The  undisturbed  ionosphere  is  a  complex  mixture  of  charged  and  neutral 
particle  species.  One  way  in  which  a  nuclear  detonation  can  modify  this  composi¬ 
tion  is  through  the  generation  of  ions  and  neutrals  with  energies  on  the  order  of 
10  or  even  larger  than  100  keV.  Once  formed,  the  fast,  high  energy  ions  can  be 
transformed  into  fast  neutrals  by  charge  exchange  collisions  while  fast  neutrals 
can  be  transformed  back  into  ions  by  impact  ionization  collisions. Between 
charge  exchange  and  impact  ionization  events,  high  energy  ions  and  neutrals  can 
increase  the  thermal  charged  particle  density  of  the  ionosphere  through  ionization 
collisions. ^  Eventually,  a  sufficiently  large  number  of  collisions  occur  so  that 
the  initially  high  energy  particles  become  part  of  the  thermal  background. 

The  persistence  of  coherent  structures,  consisting  of  fast  ions,  from 
higher  to  lower  altitudes  is  adversely  affected  by  the  charge  exchange  process. 
Ions  are  strongly  coupled  to  the  magnetic  field  geometry  but  neutrals  are  not. 
Consequently,  a  charge  exchange  collision  permits  a  fast  particle  to  traverse  a 
plasma  unimpeded  by  the  magnetic  field.  Two  different  situations  are  evident.  If 
the  probability  of  an  impact  ionization  collision  by  a  neutral  is  low  before  fast 
particles  cross  a  structural  boundary,  then  the  loss  of  fast  particles  is  convec¬ 
tive.  Conversely,  if  the  probability  is  high,  then  fast  particle  loss  is  dif¬ 
fusive.  In  general,  convective  loss  tends  to  degrade  structure  more  rapidly  than 
diffusive  loss.  Nevertheless,  if  the  neutral  mean  free  path  for  collision  is 
larger  than  the  gyroradius  of  fast  ions,  the  diffusive  loss  of  fast  particles  is 
always  larger  than  the  loss  rate  associated  with  classical  ion  diffusion  across  a 
magnetic  field. 

A  substantial  qualitative  statement,  which  is  supported  by  the  numerical 
results  subsequently  presented,  can  be  made  regarding  the  persistence  of  structure 
consisting  of  fast  ions  down  magnetic  field  lines.  Specifically,  structure 
originating  at  a  relative" v  high  alt  tude  (i.e.,  &400  km)  and  having  a  transverse 
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size  which  is  comparable  to  or  smaller  than  the  gyroradius  of  fast  ions  (e.g., 
monatomic  oxygen)  will  be  substantially  degraded  by  the  successive  charge  exchange 
and  impact  ionization  processes  before  it  reaches  lower  altitudes  (i.e., 
MOO  km).  To  demonstrate  this  assertion,  we  first  assume  that  the  structure  has  a 
sufficiently  small  size  transverse  to  the  magnetic  field  so  that  the  neutrals 
formed  following  a  charge  exchange  collision  are  convectively  lost.  This  implies 
that,  for  straight  magnetic  field  lines,  the  ion  pitch  angles  must  be  sufficiently 
oblique.  Quantitatively,  for  a  cross-section  for  impact  ionization  of  ~  4  x 
10  cm  and  a  background  neutral  density  of  —  1  x  10.  cm"  ,  the  structural  size 
perpendicular  to  the  magnetic  field  must  be  smaller  than  ~25  km.  By  comparison,  a 

3 

singly  ionized,  monatomic  oxygen  ion  with  a  speed  of  1  x  10  km/s  in  a  magnetic 
field  of  0.3  gauss  has  a  gyroradius  of  5.6  km.  When  neutrals  are  convectively  lost 
the  decay  of  structure  is  determined  by  the  rate  of  charge  exchange  collisions. 
For  this  situation,  the  structural  density  decays  in  accordance  with 


n(z)  =  n(z)  exp 


z 


(1.1) 


Here  n(z),  a,  v,  v(| ,  and  nn(z')  are  the  ion  density  at  position  z,  the  charge 
exchange  cross  section,  total  ion  speed,  ion  speed  parallel  to  the  magnetic  field, 
and  neutral  density  at  position  z',  respectively.  In  writing  Eq.  (1.1),  ax,  v,  and 
V|,  are  treated  as  constants.  For  o^lx  1Q~15  cm2,  v/vy  =1,  nn(z')  =1  x  109 
cm  ^  and  z  -  zQ  ^  100  km,  it  follows  that 


n(z)  =«  n(zQ)  exp(-lO)  =  4.5  x  10"5  n(zQ) 


(1.2) 


Equation  (1.2)  implies  severe  degradation  of  structure  with  scale  size  smaller  than 

the  fast  ion-gyroradius.  Now,  recent  fluid  simulations  of  the  EXB  gradient  drift 

instability  suggest  that  density  perturbations,  smaller  than  15  percent  on  a  larger 

1  3 

density  profile,  do  not  determine  the  outer  scale  of  the  power  spectral  density.  ’ 
It  seems  likely  that  attenuation  to  the  degree  specified  by  Eq.  (1.2)  should  be 
more  than  sufficient  to  preclude  this  structure  from  determining  the  outer  scale  at 
the  lower  altitudes. 

Prompt  striations  are  structures  with  a  scale  size  which  is  comparable  to  or 

1  4 

smaller  than  the  gyroradius  of  energetic  ions.  '  It  has  been  previously  demon¬ 
strated  that  these  structures  are  strong’/  stabilized  by  the  background  plasma  and 
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electromagnetic  effects  in  the  limit  of  small  density  gradients,  and  so  prompt 
striations  if  they  develop  are  likely  limited  to  high  altitude  (i.e.,  £.1000  km)  in 
the  HANE  environment  when  density  gradients  are  small. it  is  clear  from  the 
discussion  in  the  previous  paragraph  that  the  propagation  of  ions,  which  are 
organized  as  prompt  striations,  down  magnetic  field  lines  results  in  the  essential 
destruction  of  the  prompt  striations  once  altitudes  like  -100  to  -200  km  are 
reached. 

This  chapter  is  divided  in  the  following  way.  Section  1.2  gives  the  model 
for  the  decay  of  structure  due  to  consecutive  charge  exchange  and  impact  ionization 
collisions.  Section  1.3  presents  results  of  the  numerical  evaluation  of  the  model. 
Section  1.4  summarizes  the  results  and  conclusions. 

1.2  MODEL 

The  plasma  model  illustrated  in  Figure  1.1  is  used.  A  uniform  magnetic 
field  is  aligned  along  the  z-direction  which  is  at  an  angle  of  0, ,  with  respect  to 
gravity,  i.e., 

0J  =  cos'^g-B/dgl  |8|)]  .  (1.3) 

The  y-direction  is  ignorable  and  the  structure,  consisting  of  high  energy  ions,  has 
a  width,  2L,  in  the  x-direction.  The  calculations  assume  that  the  ions  within  the 
structure  have  the  same  speed  and  pitch  angle, 

e2  =  cos"1 ( v„ /v)  .  (1.4) 

As  the  ions  travel  down  the  magnetic  field  lines,  the  number  of  charge 
exchange  collisions  which  occur  in  an  inf initesimal  distance,  dz,  is 

dn(z)  *  n(z)Xx'1(z)  dz  ,  (1.5) 

where 

Xx(z)  ‘  2  V.1‘1 

.  n 

In  Eq.  (1.6)  the  cross  section  for  collision  with  the  background  species  is  repre 
sented  by  oxp  and  the  density  of  the  background  neutral  species  at  z  is  nn(z). 
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Now  a  fraction  of  the  charge  exchange  collisions,  f(z),  results  in  neutrals 
which  undergo  subsequent  impact  ionization  collisions  before  the  neutrals  can 
leave  the  structure.  Hence,  the  ultimate  decrease  in  the  number  of  fast  particles 
within  structure  as  a  result  of  charge  exchange  collisions  over  a  distance,  dz,  and 
a  time,  dt,  is 


dn(z)  ■  -n(z)X  -1(z)  [l  -  f(z)Jdz 


=  -n(z)Xx”i(z)  [1  -  f (z )]  V|j dt 


(1.7) 


Equation  (1.7)  implies  that 


n(z)  =  n(zQ)  exp 


/dzV‘ 


(Z')ll  -  f(z')] 


(1.8) 


Figure  1.2  illustrates  the  model  used  to  derive  an  expression  for  f(z),  the 
fraction  of  fast  neutrals  which  undergo  impact  ionization  collisions  before 
crossing  the  boundaries  of  the  structure.  The  derivation  assumes  that  the  fast  ion 
and  fast  neutral  density  within  the  structure  is  not  a  function  of  spatial  coordi¬ 
nate.  Now  charge  exchange  collisions  do  not  depend  upon  ion-gyrophase  and  so  the 
fast  neutral  velocity  distribution  function  formed  because  of  charge  exchange  is 


6(v"  -  v)  6(6  '  -  0,) 

g(v')  • - z  * 

2nv  sin(@2) 


(1.9) 


which  is  independent  of  the  azimuthal  angle  Equation  (1.9)  assumes  that  the 
fast  neutral  has  the  same  velocity  vector,  at  the  point  of  charge  exchange  col¬ 
lision,  as  the  ion  from  which  it  is  created.  It  follows  that 

Jf/2  2L  r  „  „-/£<  \ _ l 


"/C  LL 

^  J  dX  SXP  * 


x  csc(02)  csc(4>') 


(1.10) 


where 


’  2nn<z)0In 
.  n 


(1.11) 


and  «j  is  the  impact  ionization  cross  section  for  collision  of  a  fast  neutral  with 
background  neutral  species,  n.  In  writing  Eqs.  (1.10-1.11),  it  has  been  assumed 
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that  the  mean  free  path  for  impact  ionization  collision  is  evaluated  at  the  alti¬ 
tude  where  the  fast  neutrals  are  formed.  This  approximation  permits  an  analytic 
integration  of  Eq.  (1.10)  and  is  asymptotically  valid,  in  the  limit  of  greatest 
present  interest,  when  the  size  of  the  structure  transverse  to  the  magnetic  field 
is  smaller  than  the  atmospheric  scale  height.  The  integration  of  Eq.  (1.10) 
implies  that 

*  15  [al<i<a)  '  l]  +  i  *  f  a[Lo<a>Kl<a> 

*  L1(a)K0(a)]|  ,  (1.12) 

with 

a  =  [2L/Xi(z)]csc(92)  .  (1*13> 

Several  different  functions  are  used  in  Eq.  (1.12):  the  modified  Bessel  function 
of  the  second  kind  and  order  zero  [K„(a)],  the  modified  Bessel  function  of  the 
second  kind  and  order  one  [Kj(a)],  the  modified  Struve  function  of  order  zero 
[l_  (a)],*'®  and  the  modified  Struve  function  of  order  one  [L^a)].1*8 

Figure  1.3  is  a  plot  of  the  fraction,  f(z),  as  a  function  of  the  non- 
dimensional  parameter,  a.  The  figure  demonstrates  that  f(z)  has  the  correct 
asymptotic  behavior.  For  a  »1,  the  effective  mean  free  path  across  the  structure 
is  much  smaller  than  the  size  of  the  structure.  Essentially  all  the  fast  neutrals 
formed  undergo  an  impact  ionization  collision  before  they  can  cross  the  structure 
boundary  and  f(z)  -*■  1.  For  a  «  1,  the  effective  mean  free  path  across  the  struc¬ 
ture  is  much  larger  than  the  structural  size.  Fast  neutrals  are  lost  before  an 
impact  ionization  collision  can  occur  within  the  structure  and  appropriately, 
f(z)  -  0. 


1.3  RESULTS 

The  quantitative  results  presented  in  this  subsection  emphasize  monatomic 

2  3 

oxygen  particles  with  speeds  of  5  x  10  km/s  and  1  x  10  km/s.  For  a  speed  of  5  x 
102  km/s,  the  appropriate  cross  sections  for  charge  exchange  used  are  1  x  10  cm  , 
1.5  x  10“*5  cm2,  and  1.5  x  10”*5  cm2  for  collisions  of  singly  ionized,  energetic. 
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monatomic  oxygen  ions  with  monatomic  oxygen  atoms,  diatomic  oxygen  molecules,  and 

1119 

diatomic  nitrogen  molecules,  respectively.  *  ‘  *  The  cross  section  for  impact 

16  2 

ionization  is  2  x  10  cm  for  collision  of  the  energetic,  monatomic  oxygen  atoms 

with  either  monatomic  oxygen  atoms,  diatomic  oxygen  molecules,  or  diatomic  oxygen 
1119 

molecules.  ’  *  Similarly,  for  energetic,  monatomic  oxygen  ions  with  speeds  of  1 

3 

x  10  km/s,  the  cross  sections  for  charge  exchange  which  are  used  are  8  x 

162  162  162 
10  cm  ,  1  x  10"A  cm  ,  and  1  x  10  cm  ,  respectively,  for  collision  with 

monatomic  oxygen  atoms,  diatomic  oxygen  molecules,  and  diatomic  nitrogen  mole- 

1119 

cules,  respectively.  *  *  *  The  corresponding  cross  section  for  impact  ionization 
ic,  monatomic  oxygen  atoms  by  collision  with  either  monatomic  oxygen 
molecules,  or  diatomic  nitrogen  molecules  is  4  x 


of  energet 

atoms,  diatomic  oxygen 
16  9  1.1, 1.9 
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Figures  1.4-1. 7  are  plots  of  n(h)/nQ  versus  kfkm'1)  for  monatomic  oxygen  ions 
with  various  pitch  angles,  various  orientations  of  the  geomagnetic  field  with 
respect  to  gravity,  and  two  different  CIRA  reference  atmospheres.  Here  n(h)/nQ  is 
the  ion  density  at  an  altitude  h  in  units  of  the  density  at  the  starting  point  of 
the  ion  structure  which  is  assumed  to  be  at  an  altitude  of  600  km.  The  inverse  of 
the  Fourier  wavenumber,  k,  is  approximately  equivalent  to  the  structural  size 
transverse  to  the  ambient  magnetic  field.  It  is  evident  from  the  way  that  the 
figures  are  plotted  that  it  is  possible  to  determine  the  attenuation  of  structure, 
for  a  specific  wavenumber,  k,  from  some  higher  altitude,  h^,  to  some  lower  alti¬ 
tude,  by  taking  the  ratio  of  [n(h^)/nQ ]/ [n ( h2 ) /nQ ]  where  n(hj  2)/no  are  ^eter- 
mined  from  the  figures. 


Figures  1.4-1. 7  have  the  following  qualitative  form.  For  an  initial 
altitude  of  600  km,  the  disparity  in  the  attenuation  of  small  structure  as  compared 
to  large  structure  increases  as  altitude  decreases.  This  is  a  consequence  of 
particle  loss  being  more  rapid  (i.e.,  convective)  for  small  structure  and  slower 
(i.e.,  diffusive)  for  large  structure.  Because  the  convective  loss  is  only  weakly 
dependent  on  wavenumber  when  it  is  active  Icf,  Eq.  (1.1)1,  decay  is  roughly  equal 
for  all  morphology  whose  scale  size  is  smaller  than  the  mean  free  path,  across  the 
structure,  for  impact  ionization. 


Figures  1.4-1. 7  demonstrate  that  the  persistence  of  structure  down  to  lower 
altitudes  is  decreased  for  values  of  9^  deviating  from  0  and  approaching  n/2.  The 
larger  values  of  8j  imply  that  structure  originates  closer  to  the  geomagnetic 
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Figure  1. 


n(h)/n  versus  k(km"  ) /or  singly  ionized  monatomic  oxygen  ions 
with  a  speed  of  5  x  10  km/s.  The  CIRA  1972  reference  atmos¬ 
phere  for  an  exospheric  temperature  of  900°K  is  used.  The 
various  curves  from  top  to  bottom  correspond  to  altitudes  of  h  = 
500  km,  400  km,  300  km,  200  km,  and  120  km.  The  initial 
altitude  of  the  ions  is  h  =  600  km.  (Sheet  1  of  5) 
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Figure  1.5.  n(h)/n  versus  k(km"  )  /or  singly  ionized  monatomic  oxygen  ions 
with  a°speed  of  1  x  ICr  km/s.  The  CIRA  197ZL  reference  atmos¬ 
phere  for  an  exospheric  temperature  of  900°K  is  used.  The 
various  curves  from  top  to  bottom  correspond  to  altitudes  of 
500  km,  400  km,  300  km,  200  km,  and  120  km.  The  initial 
altitude  of  the  ions  is  h  =  600  km.  (Sheet  1  of  5) 


(a) 

61 

= 

n/18 

9 

92 

=  i/18 

(b) 

61 

= 

ir/18 

9 

92 

=  i/6 

(c) 

91 

= 

V18 

9 

92 

=  5 i/18 

(d) 

91 

= 

ir/18 

9 

92 

=  7i/18 

(e) 

91 

= 

ir/6 

9 

92 

=  i/18 

(f) 

91 

= 

TT/6 

9 

62 

=  i/6 

(9) 

91 

= 

TT/6 

9 

92 

=  5i/18 

(h) 

91 

s 

l/6 

9 

92 

=  7 i/18 

(i) 

91 

= 

5tt/18 

9 

92 

=  i/18 

(3) 

91 

s 

5  it/ 18 

9 

92 

*  i/6 

00 

91 

= 

5  it/18 

9 

02 

*  5i/18 

(1) 

91 

= 

5tt/18 

9 

92 

=  7 i/18 

(m) 

91 

= 

7tt/18 

9 

92 

=  i/18 

(n) 

91 

= 

7tt/18 

9 

92 

■  i/6 

(o) 

91 

= 

7  i/18 

9 

92 

•  5i/18 

(P) 

91 

= 

7  tt/18 

9 

92 

*  7i/18 

22 


Figure  1. 


n(h)/n  versus  k(km"  )  lor  singly  ionized  monatomic  oxygen  ions 
with  a°speed  of  5  x  1(T  km/s.  The  CIRA  1972  reference  atmos¬ 
phere  for  an  exospheric  temperature  of  1800°K  is  used.  The 
various  curves  from  top  to  bottom  correspond  to  altitudes  of 
500  km,  400  km,  300  km,  200  km,  and  120  km.  The  initial 


altitude 

of  the 

ions 

is 

h 

=  600  km.  (Sheet  1  of  5) 
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Figure  1.6.  Continued  (Sheet  3  of  5) 
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n(h)/n_  versus  k(km  )  for  singly  ionized  monatomic  oxygen  ions 
with  a  speed  of  1  x  10  km/s.  The  CIRA  1972  reference  atmos¬ 
phere  for  an  exospheric  temperature  of  1800°K  is  used.  The 
various  curves  from  top  to  bottom  correspond  to  altitudes  of 
500  km,  400  km,  300  km,  200  km,  and  120  km.  The  initial 
altitude  of  the  ions  is  h  =  600  km.  (Sheet  1  of  5) 
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Figure  1.7.  Continued  (Sheet  4  of  5) 
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equator.  Such  structure,  traveling  along  the  geomagnetic  field,  is  more  likely  to 
be  subject  to  a  charge  exchange  collision  as  it  goes  to  a  lower  altitude  than 
similar  structure  which  forms  at  the  geomagnetic  pole.  Stated  another  way,  between 
two  altitudes,  the  effective  columnar  density  of  neutrals  and  so  the  probability 
for  charge  exchange,  is  larger  if  the  structure  originates  nearer  to  the  geomag¬ 
netic  equator  rather  than  nearer  to  the  geomagnetic  pole. 

The  probability  that  structure  persists  to  lower  altitude  is  also  decreased 
if  the  ion  pitch  angle  approaches  */2.  Two  factors  contribute  to  this  result. 
First,  particles  with  larger  pitch  angle  must  contend  with  an  effectively  1  -ger 
columnar  density  of  neutrals  as  they  spiral  down  to  lower  altitudes.  Second, 
energetic  neutrals  formed  by  charge  exchange  are  more  likely  to  leave  the  structure 
before  they  undergo  an  impact  ionization  collision  if  the  ion  pitch  angle 
approaches  tr/2 .  This  second  factor  is  the  reason  why,  for  a^  and  &2  constants,  the 
decay  for  0^  =  a^  and  0£  =  &2  wi th  generally  results  in  larger  attenuation 
than  if  0^  =  a£  and  0£  =  aj. 

Figures  1.4-1. 5  use  a  CIRA  1972  900°K  exospheric  temperature  table  for  the 
neutral  atmosphere,  while  Figures  1.6-1. 7  use  the  1800°K  table.  The  curves  in  the 
figures  reflect  the  lower  neutral  density  which  occurs  at  a  specific  altitude  if 
the  exospheric  temperature  is  smaller.  Specifically,  structure,  especially  if  its 
transverse  size  is  smaller  than  the  effective  fast  neutral  mean  free  path,  persists 
to  a  lower  altitude  for  lower  exospheric  temperatures. 

It  is  evident  from  Figures  1.4-1. 7  that  the  combination  of  charge  exchange 
and  impact  ionization  collisions  preferential ly  destroys  fine  scale  structure 
superimposed  on  large  scale  structure.  The  result  is  a  shift  of  the  wavenumber 
spectrum  to  scale  sizes  which  are  comparable  to  or  larger  than  the  mean  free  path 
for  impact  ionization. 

Consider  a  hypothetical  structure  initially  consisting  of  energetic  ions  at 
an  altitude  above  400  km.  It  is  evident  from  the  various  curves  in  Figures  1.4-1. 7 
that  the  wavenumber  greater  than  0.1  km  components  of  this  structure  are  attenuated 
by  at  least  two  orders  of  magnitude  as  they  travel  down  to  altitudes  below  roughly 
200  km  where  the  energetic  ions  deposit  their  energy,  partially  by  ionization  of 
the  neutral  ionosphere.  At  these  altitudes,  fluid  mechanisms  like  the  EXB  gradient 
drift  instability  can  cause  fluid  structuring.  The  corresponding  outer  scale 
length  of  the  power  spectral  density  is  likely  larger  than  10  km  based  on  recent 
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fluid  simulations  of  EXB  structuring.  These  simulations  suggest  that  density 

perturbations  smaller  than  15  percent  (superimposed  on  a  larger  density  profile) 

1  3 

do  not  determine  the  long-term  outer  scale.  ’ 

The  outer  scale  length  associated  with  energetic  ionic  structure  convecting 
down  field  lines  is  at  least  10  km,  a  length  which  is  larger  than  the  expected 
gyroradii  of  monatomic  oxygen  ions.  Hence,  structures  similar  to  prompt  striations 
with  a  scale  size  comparable  to  or  smaller  than  the  ion-gyroradius*'^-*'*’  and 
originating  at  higher  altitudes  { i . e . ,  £400  km),  do  not  impose  the  long-term  outer 
scale  length  at  lower  altitudes  by  propagation  down  magnetic  field  lines.  Indeed, 
Figures  1.4-1. 7  suggest  that  structure  consisting  of  energetic  ions  may  set  an 
outer  scale  length  larger  than  100  kilometers  at  lower  altitudes,  especially  if  the 
structure  originates  near  the  magnetic  equator  and  initially  consists  of  ions  with 
large  pitch  angles. 

1.4  CONCLUDING  REMARKS 

Charge  exchange  collisions  can  convert  ions  into  neutrals,  whose  trajec¬ 
tories  are  unaffected  by  the  magnetic  field  to  lowest  order.  As  a  result,  struc¬ 
ture  consisting  initially  of  energetic  ions  can  decay  at  a  more  rapid  rate  than 
classical,  cross-field,  ion  diffusion  implicates.  This  decay  of  structure  is 
convective  or  diffusive  depending  on  whether  the  effective  neutral  mean  free  path 
for  impact  ionization  is  larger  or  smaller  than  the  structural  size  transverse  to 
the  magnetic  field.  The  degradation  of  structure  is  most  significant  when  particle 
loss  is  convective  and  particles  do  not  become  reionized  within  the  confines  of  the 
original  structural  boundaries. 

It  has  been  demonstrated  that  structures,  consisting  of  energetic  (in 
particular  monatomic  oxygen)  particles  and  having  a  scale  size  transverse  to  the 
magnetic  field  which  is  smaller  than  or  comparable  to  10  kilometers,  are  attenuated 
by  at  least  two  orders  of  magnitude  as  they  travel  down  field  lines  from  altitudes 
of  approximately  400  to  600  kilometers  down  to  altitudes  of  100  to  200  kilometers. 
Because  recent  fluid  simulations  suggest  that  density  perturbations  comparable  to 
or  less  than  15  percent,  on  a  larger  density  profile,  do  not  determine  the  long¬ 
term  outer  scale  of  the  power  spectral  density,*'3  it  seems  likely  that  the  propa¬ 
gation  of  the  structure  from  higher  to  lower  altitudes  imposes  an  outer  scale 
length  at  least  as  large  as  10  kilometers  at  the  lower  altitudes.  Origination  near 


the  magnetic  equator  and  particle  pitch  angles  approaching  90  degrees  suggest  even 
larger  outer  scale  lengths  (i.e.,  MOO  km).  Structures,  like  prompt  striations, 
which  form  at  higher  altitudes  (i.e.,  MOO  km)^"*’®  and  have  scale  sizes  perpen¬ 
dicular  to  the  magnetic  field  comparable  to  or  smaller  than  the  ion-gyroradius, 
evidently  do  not  impose  the  long-term  outer  scale  length  at  lower  altitudes  by 
propagating  down  magnetic  field  lines. 
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2.  ANALYSIS  OF  RADIATIVE  RECOMBINATION  EFFECTS  ON  DENSITY 
PROFILES  AND  THEIR  FOURIER  SPECTRUMS 


2.1  INTRODUCTION 

In  the  ionosphere,  radiative  recombination  is  an  incontrovertible  classical 

2  1-2  4 

process  which  reduces  ionization  density.  *  *  Because  it  requires  more  than  one 

particle,  a  most  significant  feature  of  the  recombination  process  is  its  inherent 

nonlinear  dependence  on  density.  Consequently,  recombination  alters  the  Fourier 

spectrum  of  density  profiles  and  so  recombination  gives  a  lower  bound  to  the 

coupling  of  Fourier  modes.  This  bound  may  be  exceeded  when  other  plasma  processes, 

like  turbulence,  are  active.  From  the  practical  standpoint,  the  change  in  Fourier 

spectrum  resulting  from  recombination  implies  an  effect  on  communications  which  is 

2  5 

in  addition  to  ionization  decay.  * 

A  simple  model  for  radiative  recombination  is  quantitatively  demonstrated 
to  cause  changes  in  both  the  morphology  and  Fourier  spectrum  of  model  one¬ 
dimensional  ionization  density  profiles.  A  key  result  of  the  analysis  is  the 
observation  that  the  short-  and  long-term  effects  of  recombination  on  structure  are 
different.  For  the  specific  density  profiles  considered,  Fourier  components  of  the 
spectrum  which  are  initially  small  have  a  general  tendency  to  at  first  increase 
relative  to  the  spatially  uniform  component.  This  effective  mode  coupling  process 
is  most  evident  for  very  cold  plasmas  (at  least  for  the  monatomic  oxygen  or  nitro¬ 
gen  model  considered  here)  in  which  the  recombination  process  is  most  nonlinear. 
Those  components  of  the  Fourier  spectrum  which  are  initially  present  decay  relative 
to  the  spatially  uniform  component.  In  the  long-term,  recombination  leads  to  the 
general  decay  of  all  Fourier  components,  resulting  from  density  inhomogeneity, 
relative  to  the  spatially  uniform  component. 

This  chapter  is  divided  into  four  sections  including  this  Introduction. 
Section  2.2  describes  the  analytic  model  used  for  the  calculations.  Section  2.3 
describes  results  from  the  numerical  evaluation  of  the  chemical  recombination 
mode.  Section  2.4  summarizes  the  major  results  and  conclusions. 
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2.2  MODEL 


As  the  basis  for  the  model,  a  weakly  ionized  neutral  plasma  is  assumed  which 

consists  of  electrons  and  singly  charged  positive  ions  recombining  into  monatomic 

atoms  of  the  same  isotopic  species.  The  evolution  of  electron  (ion)  density  as  a 

2  12  2 

function  of  time,  n(t),  is  then  determined  from  the  equation  *  * 

=  -oln(t)]  n2(t)  .  (2.1) 

In  Eq.  (2.1)  the  quantity  a(n(t)]  is  permitted  to  be  a  function  of  density.  The 
assumption, 

alnfOlsBnftj*  ,  (2.2) 


with  B  and  Y  nonnegative  constants,  permits  Eq.  (2.1)  to  be  solved  exactly  as  a 
function  of  time.  Specifically, 


V+  (Y  +  1)  Bt  n0^+1]1/(Y+ 


(2.3) 


where  nQ  is  the  initial  density  at  t  *  0. 


By  assuming  that  density  evolution  at  a  given  spatial  location  depends  on 
the  value  of  density  at  that  point  only,  it  is  evident  that  spatial  variation  can 
be  included  in  the  evolution  equation  through  the  replacement 


n(t)  -  n(x,t)  ,  nQ  >  nQ(x) 


(2.4) 


where  x  is  the  position  vector.  Equation  (2.3)  now  becomes 


n(x,t) 


nQ(?) 


[l  +  (y  *  1)  Bt  n0(J)Y+1]] 


(2.5) 


For  simplicity,  it  is  assumed  here  that  the  density  depends  on  one  spatial 
coordinate  only,  the  x-direction.  Assuming  an  arbitrary  continuous  density  vari¬ 
ation  which  is  periodic  over  a  length  2L  and  is  symmetric  about  x  e  0  permits  the 
general  density  profile  of  Eq.  (2.5)  to  be  written  in  terms  of  a  Fourier  series  in 
the  x-coordinate: 
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n(x,t)  =  —  Y/2  *  S  ~  P'"l/2  cos  ^  x^ 
(2L)1/2  ~  (L)1/Z  L 


(2.6) 


In  Eq.  (2.6)  the  time  dependent  Fourier  amplitudes  are  determined  by 

L 

aft)  =  - f - t-pt — T - 7 - rr t  f  dx  x)  »  (2-7) 

with  6  the  Kronecker  delta.  It  follows  from  Eq.  (2.7)  that  the  quantity, 
1 

(21)  7  aQ(t),  is  the  total  ion  density  between  x  =  -L  and  x  =  +L  at  any  time  t  and 
is  the  only  Fourier  amplitude  present  in  an  uniform  plasma. 

It  is  useful  to  examine  the  form  of  the  Fourier  amplitudes  in  two  limits 
[i.e.,  pf+1)  6t  nQ(x)Y+^«  1  and  »lj.  For  (Y+l)  8t  n^x^^cd  for  all  x,  the 
binomial  expansion  permits  Eq.  (2.7)  to  be  approximately  written 


aP(t’  %  ll  ♦  [<2)I/2-l]Vo  I  WV2 


L 

/dx  n 

0 


(x) 


|l  -  fit  n0(x)Y+1|  cos  (p  x) 


Y+1 

Similarly,  for  (Y+l)  8t  no(x)  A»  1  for  all  x. 


a„(t) 


P  +  [(2)1/2-l]<5p>oi  (L) 


i  1/2 


L 

/ 


dx 


cos(pirx/L) 

TTTf+TJ 


I(Y+1)  fit] 


1--^! 

Bt  n0(x)Y+1 


(2.8) 


(2.9) 


Equation  (2.9)  implies  that. 


«n(t)  * 


(2L) 


1/2 


’(Y+1)  fit]17^*^ 


and 


p*0 


_ 2  l/(Y+l)(2T*3)/(Y4l)  d  cqs(pttx/L) 

< (8t)<Y+^Y+1>  J  nJx)T+1 


(2.10a) 


(2.10b) 
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Explicit  quantitative  temporal  and  spatial  evaluation  of  the  density  pro¬ 
file  when  recombination  is  active  requires  specific  values  for  the  parameters,  8 
and  Y.  In  this  paper  it  is  assumed  that  the  ion  and  neutral  species  are  monatomic 
nitrogen  or  oxygen.  The  particular  exothermic  reactions  considered  are: 


N+  +  e”  -*■  N  +  14.53  eV 
0+  +  e"  -*•  0  +  13.61  eV 


(2.11) 


where  N+,  e  ,  N,  0+,  and  0  represent  a  positive  nitrogen  ion,  an  electron,  a 

neutral  nitrogen  particle,  a  positive  oxygen  ion,  and  a  neutral  oxygen  particle, 

respectively.  For  the  reactions  specified  by  Eq.  (2.11),  the  values  of  a(n)  in 

Eq.  (2.1)  are  tabulated  in  Table  2.1  for  different  densities  and  for  electron 

?  6 

temperatures,  T  ,  of  0.043  eV  and  0.518  eV,  respectively.  '  A  linear  regression 
analysis  for  Hn[a(n)J  gives 

Te  =  0.043  eV  ,  8  = 4.433  x  10“17  ,  Y  ^  0.6689  (2.12a) 

Te  =  0.518  eV  ,  8 35  1.202  x  10‘13  ,  Y^0. 04565  (2.12b) 


with  8  in  cgs  units.  Clearly,  the  recombination  coefficient  ot(n)  is  more  sensi¬ 
tively  dependent  on  density  for  lower  temperatures  than  for  higher  temperatures. 


2.3  EXAMPLES 

To  demonstrate  the  quantitative  decay  of  ionized  density  structure  as  a 
function  of  time,  four  different  initial  density  profiles  are  considered: 


nQ(x)  =  N [l  -  0.3  (x/l)]  ,  (2.13a) 

•  N[l  -  0.7  (x/L)]  ,  (2.13b) 

*  N  |l  +  0.1  [cos(4jtx/L)  +  cos(5ttx/L)]{  ,  (2.13c) 

=  Njl  +  0.3  [cos(4irx/L)  +  cos(5ttx/L)]  (  .  (2.13d) 
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Table  2.1.  Values  of  a(n)  in  units  of  cm  /s  for  various  densities  and 
electron  temperatures,  T  ,  of  0.043  eV  and  0.518  eV 


Te(eV) 
n(cnf  3) 

0.043 

0.518 

107 

4.1  x  10‘12 

2.5  x  10'13 

108 

9.3  x  ID'12 

2.7  x  10'13 

109 

2.7  x  10'11 

3.1  x  10'13 

o 

rH 

o 

1.1  x  10"10 

3.6  x  10'13 

1011 

1.0  x  IQ'9 

4.0  x  10‘13 

1012 

9.0  x  10’9 

4.0  x  10'13 

In  Eqs.  (2.13a)-(2.13d)  the  parameter  N  is  a  constant  which  is  determined  by  the 
actual  quantitative  density  at  some  specified  point  in  the  density  profile  (e.g.,  x 
*  0)  and  2L  is  the  periodicity  length.  Figures  2.1a-d  are  plots,  for  the  four 
different  density  profiles,  of  a0(t)/N  versus 

e «  n[(Y+1)  &t]1/^Y+1^  ,  (2.14) 


for  recombination  at  Tg  =  0.043  eV  and  Tg  =  0.518  eV.  Because  aQ(t)  is  directly 
proportional  to  the  total  ion  density  for  -L  <  x  <  L  at  any  time,  t,  it  is  evident 
from  the  figures  that  the  rate  of  total  density  decay  for  the  two  values  of  Tg 
considered,  is  most  rapid  for  e  <  1  when  density  is  greatest.  This  feature  of  the 
figures  is  clearly  consistent  with  Eq.  (2.1)  and  the  positive  values  of  B  and  Y 
given  by  Eqs.  (2.12a-b).  There  are  two  other  ubiquitous  features  of  Fig¬ 
ures  2.1a-d.  First,  the  curves  for  Tg  =  0.043  eV  decrease  less  rapidly  than  those 
for  Tg  =  0.518  eV  when  e  is  small.  This  characteristic  of  the  curves  is,  from 
Eqs.  (2.8)  and  (2.14),  a  consequence  of 

aQ(t)  *  (2L)1/2  N  J  dx 

o 

0 

and  Y  being  larger  for  Tg  =  0.043  eV  than  for  Tg  *  0.518  eV.  Second,  for  asymp¬ 
totically  large  values  of  e  the  curves,  for  the  two  different  temperatures, 
approach  each  other  because 


1 

Y+l 


m 


(2.15a) 


,(t>  „ 


(2.16a) 


da0(t)  (2L)1/2  N 
de  '  "  -2 


(2.16b) 


which  is  independent  of  Y. 

Figures  2.2a-d  and  2.3a-d  are  plots  of  the  four  density  profiles  in  normal¬ 
ized  units,  n(x,t)/N,  as  a  function  of  (x/L)  for  Tg  *  0.043  eV  and  0.518  eV, 
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Figure  2.1a.  a0(t)/N  versus  e  for  Te  =  0.043  eV  and  0.518  eV, 

respectively.  Initial  density  profile  specified 
by  Eq.  (2.13a). 
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1.0  2.0  3.0  4.0  5.0  6.0  7.0 


Figure  2.1b.  a0(t)/N  versus  e  for  T  =  0.043  eV  and  0.518  eV, 
respectively.  Initiardensity  profile  specified 
by  Eq.  (2.13b). 


47 


a0(t)/N  versus  e  for  Te  =  0.043  eV  and  0.518  eV, 
respectively.  Initial  density  profile  specified 
by  Eq.  (2.13c). 


aQ(t)/N 


1.0  2.0  3.0  4.0  5.0  6.0  7 


Figure  2. Id.  a0(t)/N  versus  e  for  Te  =  0.043  eV  and  0.518  eV, 
respectively.  Initial  density  profile  specified 
by  Eq.  (2.13d). 
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(x)/N  =1-0.3  (x/L) 


profile  specified  by  Eq. 


(x)/N  =1-0.7  (x/L) 


Figure  2.2b.  n(x,t)/N  versus  x/L  for  Tp  =  0.043  eV  and  various 


x)/N  -  1.0  +  0.1  [co$(4hx/L)  +  cos(5nx/L)] 


for  Te  =  0.043  e 
ile  specified  by 


(x)/N  =1-0.3  (x/L) 


518  eV  and 
ed  by  Eq. 


(x)/N  =1-0.7  (x/L) 


n(x,t)/N  versus  x/L  for  Te  =  0.518  eV  and 
Initial  density  profile  specified  by  Eq. 


profile  specified  by 


n(x,t) 


x/L  for  Tg  =  0.518  eV  and  various 
profile  specified  by  Eq.  (2.13d). 


respectively.  In  the  figures  the  various  curves  correspond  to  different  assumed 
values  of  e.  A  general  characteristic  of  all  the  figures  is  the  tendency  for 
recombination  to  simultaneously  reduce  both  the  absolute  density  and  density 
inhomogeneity.  This  characteristic  is  a  consequence  of  the  more  rapid  decay  of 
charged  particle  density  in  spatial  locations  where  the  density  is  initially 
largest  and  the  weaker  decay  at  locations  where  density  is  initially  smallest. 

Because  the  recombination  decay  rate  is  strongly  dependent  upon  density, 
the  Fourier  spectrum  of  the  various  initially  inhomogeneous  density  profiles  must 
be  a  function  of  time.  Indeed,  Figures  2.4a-d  and  2.5a-d,  which  are  plots  of 
lap(t)/a0(t)l  versus  e  for  different  values  of  harmonic  number,  p,  clearly  demon¬ 
strate  the  time  variation  of  the  Fourier  coefficients  for  the  various  density 
profiles  considered.  In  Figures  2.4a-d  and  2.5a-d  the  temperatures  considered  are 
Tg  =  0.043  eV  and  0.518  eV,  respective’/-  The  figures  have  several  significant 
features.  First,  the  absolute  magnitude  of  |ap(t)/aQ(t)|  is  dependent  upon  the 
degree  of  density  inhomogeneity.  Specifically,  as  inhomogeneity  increases  the 
contribution  of  the  p  f  0  Fourier  amplitudes  which  .  -e  initially  present  always 
decay  relative  to  aQ(t).  This  decay  is  a  consequence  of  the  density  inhomogeneity 
erosion  associated  with  recombination.  Third,  the  p^O  Fourier  amplitudes  which 
are  not  initially  present  generally  grow  relative  to  aQ(t)  for  small  values  of  e. 
Only  Figure  2.4d  shows  evidence  of  more  complicated  behavior  superimposed  on  the 
general  growth.  The  relative  growth  is  caused  by  both  the  change  in  density 
profile  shape  and  the  simultaneous  general  decay  of  aQ(t)  ( i . e . ,  total  density). 
The  more  complicated  feature  of  Figure  2.4d  is  a  consequence  of  the  change  in  sign 
of  certain  ap(t)  as  e  increases  from  zero.  Fourth,  those  |ap(t)/aQ(t)|  which 
initially  increase  attain  their  maximum  value  for  e  ~  1.  For  larger  c,  these 
|ap(t)/a0(t)|  decrease.  Fifth,  for  asymptotically  large  values  of  e  al  1  la  (t)/ 
aQ(t)|  for  a  specific  electron  temperature  decay  at  the  same  rate.  In  particular 
from  Eqs.  (2.10a-b), 

l*p(t)/a0(t)l  ~  l/eY+1  .  (2-17) 

is  independent  of  density  profile  and  so  density  profiles  become  uniform  on  a 
faster  time  scale  than  total  ionization  decay.  Because  the  value  of  Y  is  larger 
for  Tg  *  0.043  eV  than  for  Tg  *  0.518  eV,  the  relative  decrease  of  the  Fourier 
amplitudes  as  specified  by  Eq.  (2.17)  is  more  rapid  at  the  lower  temperature. 
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1.0  2.0  3.0  4.0  5.0  6.0 


Figure  2.4a. 


ap(t)/a0(t)|  versus  e  for  T-  =  0.043  eV  and  various  p. 
nitial  density  profile  specified  by  Eq.  (2.13a) 
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ap(t)/a0(t)i 


Figure  2.4b.  |ap(t)/aq(t) |  versus  e  for  T*  =  0.043  eV  and  various 
Initial  density  profile  specified  by  Eq.  (2.13b). 
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j(t)/a0(t)|  versus  e  for  T_  =  0.043  eV  and  various  p 
itial  density  profile  specified  by  Eq.  (2.13c). 


Figure  2.4d.  |ap(t)/a0(t) |  versus  e  for  Tg  =  0.043  eV  and  various  p 

Initial  density  profile  specified  by  Eq.  (2.13d). 


Vt)/ao(t)| 


c 


Figure  2.5a.  |aD(t)/ag(t) j  versus  e  for  Tg  =  0.518  eV  and  various  p. 

Initial  density  profile  specified  by  Eq.  (2.13a). 
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|an(t)/an(t) 


Figure  2.5c. 


|  versus  z  for  Te  =  0.518  eV  and  various  p 
ity  profile  specified  by  Eq.  (2.13c). 


Figure  2.5d.  |ap(t)/a0(t) j  versus  c  for  T?  =  0.518  eV  and  various  p 

Initial  density  profile  specified  by  Eq.  (2.13d). 
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I 

t 

j 

The  qualitative  variation  of  certain  of  the  |ap ( t ) /a0( t ) |  on  e  has  been  shown 

|  to  be  different  depending  on  whether  e  <1  or  3,1.  This  change  in  behavior  can  be 

j  related  to  real  time  by  using  Eqs.  (2.10a-b)  and  (2.14).  Figures  2.6a-b  are  plots 

of  t(second)  versus  e  for  T  =  0.043  and  0.518  eV,  respectively.  The  different 
1  e  -3 

|  curves  in  each  of  the  figures  correspond  to  different  choices  of  density,  N(cm  ). 

The  curves  readily  show  that  the  recombination  time  scale,  defined  as  the  time  when 

e=  1,  is  smallest  when  density  is  large  and  temperature  is  low.  For  example,  in 

9  -3 

i  Figure  2.6a,  e  =  1  at  approximately  12  seconds  for  N  =  1  x  10  cm  and  at  approxi- 

i  n  7-3 

mately  2.8  x  10  seconds  for  N  =  1  x  10  cm  .  Similarly  in  Figure  2.6b,  e  =  1  at 

11-3  3 

approximately  25  seconds  for  N  =  1  x  10  cm  and  at  approximately  3  x  10  seconds 
g  >3 

for  N  =  1  x  10  cm  .  Other  values  of  e  and  the  times  to  which  they  correspond  can 
|  be  read  from  the  figures. 

2.4  SUMMARY  AND  CONCLUSIONS 

Recombination  is  a  classical  process  which  can  both  reduce  total  ionization 
and  smooth  out  ionization  density  profiles  irrespective  of  the  plasma  turbulence 
level.  Hence,  recombination  is  one  classical  process  which  provides  a  lower  bound 
on  the  level  of  turbulence  necessary  to  make  the  plasma  behavior  anomalous. 

Because  the  temperature  of  the  ionosphere  is  often  comparable  to  0.1  eV  and 

the  primary  ionic  species  in  the  ionosphere  is  monatomic  oxygen,  the  radiative 

2  1-2  4 

combination  is  a  relevant  process  for  consideration  in  the  ionosphere.  ' 

!  Model  calculations  of  the  radiative  recombination  process  for  several  different 

i  one-dimensional  density  profiles  have  been  described  for  a  neutral  plasma  con¬ 

taining  monatomic  oxygen  atoms,  the  main  neutral  constituent  in  the  ionosphere 

2  1 

above  approximately  300  kilometers.  ’  Qualitatively  and  quantitatively  the 

2  6 

radiative  recombination  process  is  identical  for  monatomic  nitrogen. 

The  evaluation  of  the  nonlinear  recombination  equation,  Eq.  (2.1),  demon¬ 
strates  the  simultaneous  change  in  the  shape  and  total  ionization  of  the  density 
profiles.  Time  variation  in  the  relative  size  of  Fourier  amplitudes  is  a  direct 
j  consequence  of  the  change  in  density  morphology.  In  particular,  for  very  short 

times,  such  that  e«l,  those  Fourier  amplitudes  not  initially  present  for  the 
density  profile  have  a  general  tendency  to  grow  relative  to  aQ(t).  Subsequently, 
for  e»l,  all  p  f  0  Fourier  amplitudes  decay  in  the  same  functional  way  relative  to 
aQ(t).  For  long-times,  aQ(t)  and  total  ionized  particle  density  decay  inversely 
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with  time  to  a  power  which  depends  on  the  electron  temperature.  The  indicated 
functional  behavior  of  long-term  Fourier  amplitude  decay  with  time  is  a  signature 
of  the  radiative  recombination  process. 
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3.  NONADIABATI CITY  AND  THE  HIGH-ALTITUDE  NUCLEAR  ENVIRONMENT 


3.1  INTRODUCTION 

The  magnetic  moment  is  an  adiabatic  constant  of  motion  in  magnetic  field 

configurations  with  sufficiently  gentle  magnetic  field  gradients  and  slow  temporal 

variations.  However,  if  spatial  inhomogeneities  in  the  magnetic  field  are  large 

enough  or  if  the  magnetic  fields  change  fast  enough,  the  magnetic  moment  is  not  a 

reasonable  constant  of  the  motion.  In  fact,  numerical  studies  of  particle 

motion  demonstrate  a  rather  sharp  transition  between  adiabatic  and  nonadiabatic 

3  1-3  4 

behavior  as  the  parameter  X  =  P/L  is  increased.  *  ’  Here,  P  =  w/Q.  is  the 

particle  gyroradius  for  speed,  v,  and  gyrofrequency,  Q,  while  L  is  the  inhomo¬ 
geneity  scale  length.  Particles,  whose  motion  does  not  conserve  magnetic  moment, 
generally  have  nonadiabatic  changes  in  their  magnetic  moment  as  they  pass  through 
points  of  minimum  (maximum)  magnetic  field  strength  on  their  guiding  center 
trajectories.  The  magnitude  of  these  nonadiabatic  jumps  depends  on  the  particle 
pitch  angle  and  gyration  phase. 

There  are  two  situations  where  the  nonadiabatic  particle  motion  can  in¬ 
fluence  the  ionospheric  and  magnetospheric  environments  following  a  high-altitude 
nuclear  burst.  First,  a  nuclear  burst  can  push  aside  the  ambient  geomagnetic  field 

and  temporarily  create  a  new  magnetic  field  configuration  containing  local 

3  5  3  6 

extremums  in  magnetic  field  strength.  ’  '  Superthermal  ions  flowing  down  the 

field  lines  in  this  configuration  can  undergo  nonadiabatic  jumps  in  magnetic  moment 

as  they  traverse  these  extremums.  As  a  result,  initial  structure  whose  size  is 

smaller  than  the  maximum  ion  gyroradius  is  fuzzed  on  a  size  scale  comparable  to 

this  gyroradius.  Hence,  nonadiabaticity,  like  charge  exchange,  supports  the 

premise  of  outer  scale  length  in  the  power  spectral  density  being  at  least  as  large 

3  6 

as  the  gyroradii  of  fast  ions.  Second,  nuclear  bursts  may  inject  large  quanti- 

3  7 

ties  of  superthermal  ion  material  into  the  magnetosphere.  However,  only  ions, 
whose  trajectories  sensibly  conserve  magnetic  moment,  can  be  confined  on  drift 
orbits  in  the  magnetosphere  and  contribute  to  a  magnetospheric  ring  current.  These 


ring  currents  not  only  have  sharp  gradients,  which  may  directly  influence  electro¬ 
magnetic  transmissions  from  satellites,  but  by  their  inherent  nature  can  affect  the 
whole  ionospheric  and  magnetospheric  current  system.  ’  For  example,  associated 
field-aligned  currents  may  drive  conventional  ionospheric  structuring  mechanisms 
like  the  current-convective  instability.  Also,  electron  currents  may  permit 
shorting  of  ambipolar  electric  fields  in  the  ionosphere  and  a  concomitant  increase 
in  the  basic  ionospheric  plasma  diffusion  rate.  Changes  in  the  rate  of  plasma 
diffusion  can  affect  striation  morphology. 

This  chapter  shows  the  relevance  of  nonadiabatic  scattering  of  ions  to  the 
phenomenology  of  the  high-altitude  nuclear  environment.  Section  3.2  gives  a  semi- 
quantitative  description  of  conditions  required  for  nonadiabatic  changes  in 
magnetic  moment.  Also,  a  Fokker-Planck  model  is  described  which  is  appropriate  to 
situations  where  the  jumps  in  magnetic  moment  are  relatively  small  compared  to  the 
initial  magnetic  moment.  Section  3.3  shows  the  results  of  an  exact  numerical 
evaluation  of  fast  ion  trajectories  for  a  static  magnetic  field  configuration  with 
gradients  which  can  be  comparable  to  those  associated  with  a  high-altitude  nuclear 
burst.  These  calculations  clearly  demonstrate  that  nonadiabatic  jumps  in  magnetic 
moment  can  occur  in  an  appropriately  inhomogeneous  magnetic  field  configuration. 
Section  3.4  shows  that  superthermal  ions  generated  by  a  high-altitude  nuclear 
burst  can  undergo  nonadiabatic  changes  in  magnetic  moment  in  the  magnetosphere. 
Section  3.5  summarizes  the  results  and  gives  concluding  remarks. 

3.2  FOKKER-PLANCK  MODEL  FOR  NONADIABATICITY 

Numerical  studies  of  particle  motion  in  complex  magnetic  geometries  show 

3  1-3  4 

that  particle  trajectories  begin  to  evince  nonadiabatic  behavior  when  ' 

v/n  >  0.04  L  ,  (3.1) 

3  2  3  4 

and  demonstrate  large,  nonperturbative  changes  in  magnetic  moment  when  '  ’  ’ 

v/fl  £  0.1  L  .  (3.2) 

O 

For  singly  ionized  monatomic  oxygen  ions  with  a  speed,  v  s  1  x  10  cm/s,  and  a 
magnetic  field  strength  of  0.3  gauss,  Eqs.  (3.1)  and  (3.2)  imply  that  L  5. 140  km 
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for  the  onset  of  nonadiabaticity  and  L  5  56  km  for  very  large  changes  in  magnetic 
moment. 

Whereas  the  magnetic  moment,  u,  can  undergo  jumps  as  particles  traverse 
extremums  in  magnetic  field  strength,  the  sign  of  the  jump,  Au,  depends  on  the 
gyrophase  of  the  electron  at  the  minimum-  (maximum)  point.  Hence, 


<Ap>  *  0 


(3.3) 


where  the  angle  bracket  denotes  an  average  over  gyrophase.  The  mean  square 
quantity,  <Ap>,^'2  averaged  over  gyrophase,  can  be  expressed  as3«l-3.4 


vlo 


(3.4) 


where  viQ  is  the  perpendicular  particle  speed  at  the  minimum  (maximum)  point  of 

magnetic  field  strength.  The  coefficient  A  can  be  taken  to  be  roughly  a  constant, 

4 .3. 1-3. 4  jn  generai ^  K(v2/viQ2)  has  a  detailed  complicated  functional  dependence 

which  depends  on  geometry  but  the  qualitative  features  of  the  function  has  the  form 

3  13  4  7  7 

illustrated  in  Figure  3.1.  *  The  value  of  K(v  /v  )  is  a  strong  function  of 

partide  pitch  angle,  especially  for  trapped  particles,  v  /v  <  2.  The  value  of 
K(v2/vln2)  increases  with  the  pitch  angle,  i.e.,  K(v2/v,  2)  »1  for  (v/v,J2  -M. 
K(v  /vi0  )  is  roughly  constant  (of  the  order  of  1)  for  the  regime  (v/vi  )4  &4. 
Since  for  local  perpendicular  speed,  v^  local  magnetic  field  strength,  B,  and 
particle  mass,  m. 


mv, 


u  = 


~W 


(3.5) 


the  change  in  pitch  angle  described  by  Eq.  (3.4)  goes  to  zero  algebraically  for 

small  pitch  angle  and  exponentially  for  large  pitch  angle.  Nonadiabatic  jumps  in 

3  13  4 

magnetic  moment  require  oblique  pitch  angle.  ’  *  ’ 

A  velocity  space  diffusion  operator  has  previously  been  developed  for  the 

3  7 

nonrelativistic  limit  using  a  Fokker-Planck  formalism,  appropriate  to  the  limit 
2  2 

<(Au)  >/u  «  1.  According  to  Reference  3.2,  the  diffusion  operator  describing 
nonadiabaticity  effects,  N,  has  the  following  form. 
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(3.6) 


with  f,  the  velocity  distribution  function.  For  trapped  particles,  At  is  the 
bounce  period,  while  for  untrapped  electrons.  At  is  the  time  required  to  traverse 
the  magnetic  field  geometry  of  interest. 


For  application  to  the  high-altitude  nuclear  environment,  it  is  useful  to 
rewrite  Eq.  (3.6)  in  terms  of  particle  momentum  components  at  the  points  which 
correspond  to  extremums  in  magnetic  field  strength.  In  particular. 


with 


N(Pi0’pJf 


3  3  \  <(Au)2>  "2Bo2 

3pi0  3PiJ  2At  ^7" 
3P^  "  3P^)  f(Pio’P|,o) 


(3.7) 


Pio  3  mv; o  •  Pjio  *  "^tio  *  (3‘8) 

and  the  parallel  velocity,  v||0,  at  the  point  with  the  minimum  (maximum)  magnetic 
field  strength,  BQ. 

An  H-theorem  can  be  derived  from  Eq.  (3.7): 

°->/ d3P  f(Pa0’P|io)  N{Pio'Piio)  f^io*p-.o) 


/ 


d3P 


<(AU)^> 

2At 


m2B  2 
o 

~~T 


r/pHo  3  _  3  \ 

IVlO  3P10  3P|io/ 


,  2 


f(Pio’Piio> 


(3.9) 


From  Eq.  (3.9)  it  follows  that  for  the  general  case  of  initially  oblique  particle 
pitch  angles,  the  nonadiabatic  jumps  in  magnetic  moment  generally  tend  to  drive  the 
distribution  function  to  Isotropy  while  conserving  energy.  There  are  also  special 
steady-state  distribution  functions  corresponding  to  zero  and  ninety  degree  pitch 
angles,  respectively. 


75 


In  plasmas  with  rotational  or  translational  symmetry  there  is  a  corres- 

3  8 

ponding  canonical  momentum  which  is  a  constant  of  motion.  This  canonical 
momentum  determines  the  trajectory  in  real  space  which  can  be  traversed  by  the 
gyrocenter  of  a  specific  particle.  Nonadiabatic  changes  in  magnetic  moment  must 
preserve  constant  canonical  momentum  and  so  nonadiabaticity,  by  itself,  does  not 
necessarily  imply  particle  diffusion  across  magnetic  fields.  However,  nonadia¬ 
baticity  can  still  result  in  the  loss  of  particles  by  pushing  them  into  loss  cones. 
For  systems  which  do  not  have  a  coordinate  of  symmetry,  the  kinetic  diffusion 
operator,  Eq.  (3.7),  has  to  be  written  in  a  more  general  form  in  order  to  accu¬ 
rately  assess  the  diffusion  associated  with  nonadiabatic  changes  in  the  magnetic 
moment. 


3.3  EXAMPLE  FOR  THE  VICINITY  OF  AN  IONOSPHERIC  NUCLEAR  BURST 

In  this  section  we  present  numerical  results  which  demonstrate  that  non¬ 
adiabaticity  is  a  relevant  process  for  magnetic  field  geometries  with  gradients 
comparable  to  those  which  might  be  expected  in  the  vicinity  of  the  magnetic  bubble 
formed  following  an  ionospheric  nuclear  burst.  For  these  examples  it  was  numer¬ 
ically  convenient  to  use  a  simple  mirror  geometry  to  demonstrate  the  nonadiabatic 

changes  in  magnetic  moment.  The  particular  geometry  chosen  for  the  present 

3  9 

analysis  is  appropriate  to  a  straight,  vacuum,  mirror  field:  ’ 


B  (r,z)  =  B  [a-  IQ(p.  r)  cos(^-  z)] 

o  Lo 

(3.10a) 

B  (r,z)  =  -  B0  IjfA  r)  sin(4-  z) 

L0  0 

(3.10b) 

Here  a,  LQ,  and  BQ  are  a  numerical  parameter,  a  scale  length,  and  a  magnetic  field 
strength,  respectively.  These  quantities  are  chosen  to  have  values  consistent  with 
the  high-altitude  nuclear  environment.  The  symbols,  r  and  z,  are  the  radial  and 
axial  coordinates  of  the  straight  cylinder,  respectively.  The  functions,  IQ  and 
Ip  are  the  modified  Bessel  functions  of  the  first  kind  and  of  order  zero  and  one, 
respectively.  The  form  for  the  fields  assumes  that  only  vacuum  magnetic  fields  are 
present  and  neglects  plasma  diamagnetism. 
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Figures  3.2a-b,  which  are  plots  of  2y/mc  (G  )  versus  z(km),  are  a  conse¬ 
quence  of  the  solution  to  the  exact  equation  of  motion  in  the  magnetic  field 
geometry  specified  by  Eqs.  (3.10a-b).  The  speed  of  light  is  represented  by  c.  The 
values  of  ot(13),  LQ  (20  km),  and  BQ  (0.027  G)  actually  imply  minimum  and  maximum 
magnetic  field  strengths  on  axis  (i.e.,  r  =  0)  of  -0.32  and  ~0.38  G,  respec¬ 
tively.  The  particle  considered  is  singly  ionized  monatomic  oxygen  with  an  initial 
pitch  angle  of  45°  and  an  initial  location  of  r  =  10  km  and  z  =  0  km.  The  magnetic 
field  strength  at  the  initial  particle  location  is  0.3  G. 

3 

Figure  3.2a,  which  assumes  an  ion  speed  of  1  x  10  km/s,  shows  that  the 
magnetic  moment  is  not  an  adiabatic  constant  for  the  specified  parameters.  The 
initial  ion  position  is  at  z  =  -20  km  and  the  initial  ion  velocity  is  in  the 
positive  z-direction.  Away  from  z  =  0,  there  are  oscillations  in  the  magnetic 
moment  which  are  directly  related  to  the  ion  gyration  about  the  magnetic  field 
lines.  Near  the  point  of  minimum  magnetic  field  strength  at  z  =  0,  there  is  a 
nonadiabatic  increase  in  magnetic  moment.  Just  before  z  =  20  km,  the  particle's 
parallel  velocity  becomes  zero  and  undergoes  a  single  bounce  motion  (not  shown  in 
Figure  3.2a)  before  crossing  the  z  =  20  km  point.  It  is  evident  from  Figure  3.2a 
that  a  particle  which  is  initially  untrapped  can  become  temporarily  trapped  around 
a  point  of  minimum  magnetic  field  strength  if  the  magnetic  moment  is  not  constant. 
Conversely,  a  particle  which  is  initially  trapped  can  become  untrapped  if  the 
magnetic  moment  is  not  an  adiabatic  invariant.  Such  behavior  is  impossible  for 
situations  where  the  magnetic  moment  is  sensibly  constant. 

Figure  3.2b  assumes  the  same  parameters  as  Figure  3.2a  except  that  the  ion 
speed  is  2  x  10^  km/s.  The  particle  motion  is  initially  in  the  positive  direction 
and  begins  with  the  solid  line  at  z  *  -20  km.  As  in  Figure  3.2a,  there  are 
oscillations  in  the  magnetic  moment  away  from  z  =  0.  Near  z  =  0  there  is  a  mild 
nonadiabatic  change  in  the  magnetic  moment  but  the  ion  continues  to  travel  to  z  = 
20  km.  The  particle  trajectory  corresponding  to  the  dashed  line  begins  at  z  = 
-20  km  with  initial  ion  velocity  components  identical  to  the  final  ones  at  z  = 
20  km  for  the  solid  line.  The  dashed  line  demonstrates  a  dramatically  different 
behavior  for  the  magnetic  moment  than  the  solid  line.  This  line  shows  a  substan¬ 
tial  nonadiabatic  increase  in  magnetic  moment  near  the  point  of  minimum  magnetic 
field  strength  at  z  *  0.  Beyond  z  *  0,  the  ion  has  a  subsequent  reversal  in 
trajectory  and  the  ion  moves  in  the  negative  direction.  Only  part  of  this  motion 
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in  the  negative  direction  is  shown  in  Figure  3.2b.  The  trapping  of  the  ion  by  the 
magnetic  field  is  only  temporary  and  eventually  the  ion  crosses  the  2  =  20  km  point 
after  one  bounce  motion.  In  addition  to  demonstrating  the  possibility  of  both 
trapped  and  untrapped  particle  dynamics  when  the  magnetic  moment  is  not  conserved. 
Figure  3.2b  shows  that  the  nonadiabatic  change  in  magnetic  moment  is  a  function  of 
pitch  angle  and  gyrophase.  This  result  follows  from  the  different  trajectories 
associated  with  the  solid  and  dashed  lines. 

During  early  times,  nuclear  bursts  in  the  ionosphere  strongly  modify  the 

o  c  ^  c 

ambient  magnetic  field  geometry,  for  example,  by  creating  magnetic  bubbles.  * 
Energetic  ions  traveling  through  this  geometry  may  cross  points  of  minimum  or 
maximum  magnetic  field  strength  with  resulting  substantial  changes  in  magnetic 
moment  and  pitch  angle.  The  nonadiabatic  jumps  in  magnetic  moment  are  a  function 
of  pitch  angle  and  gyrophase. 

Consider  a  structure  consisting  of  energetic  ions  and  having  a  size  trans¬ 
verse  to  the  magnetic  field  which  is  initially  smaller  than  the  maximum  gyroradius 
possible  for  a  given  speed.  It  is  evident  that  nonadiabaticity  tends  to  broaden 
the  transverse  size  and  make  it  comparable  to  the  maximum  gyroradius.  Nonadia¬ 
baticity  supports  the  premise  of  minimum  outer  scale  length  being  comparable  to  the 

3  6 

early-time  energetic  ion-gyroradius. 

3.4  NONADIABATICITY  IN  THE  MAGNETOSPHERE 

Following  a  high-altitude  nuclear  burst,  ions  may  be  injected  into  the 
magnetosphere.  Only  those  ions,  with  reasonably  constant  magnetic  moments,  can  be 
confined  on  drift  orbits  in  the  magnetosphere  and  contribute  in  a  substantial  way 
to  the  magnetospheric  ring  current.  These  currents,  by  virtue  of  their  sharp 
gradients,  may  not  only  directly  affect  transmissions  from  satellites  but  also  may 
have  a  direct  impact  on  the  whole  ionospheric  and  magnetospheric  current  system. ^ 
Associated  field-aligned  currents  may  drive  conventional  ionospheric  structuring 
mechanisms  like  the  current-convective  instability.  Also,  electron  currents  may 
permit  shorting  of  ambipolar  electric  fields  in  the  ionosphere  and  permit  a 
concomitant  increase  in  the  basic  ionospheric  plasma  diffusion  rate  and  the  inner 
scale  length  of  the  power  spectral  density. 

The  consideration  of  nonadiabatic  particle  behavior  permits  a  bound  to  be 
put  on  the  maximum  energy  of  ions  confined  to  the  magnetosphere.  This  bound  does 
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not  depend  on  a  hypothetical  anomalous  mechanism  but  rather  only  depends  on  the 
inherent  magnetic  field  geometry  of  the  magnetosphere. 

To  demonstrate  that  nonadiabaticity  is  a  valid  consideration  for  ions 
injected  into  the  magnetosphere  following  a  nuclear  burst,  we  consider  a  dipole 
magnetic  field  geometry  for  the  earth^'^ 


BR 


2M  sin(em) 


(3.11a) 


M  cos(e  ) 


(3.11b) 


■} 

with  M  =  8  x  10  G-cm  .  The  radial  distance  from  the  center  of  the  earth  and  the 
magnetic  latitude  are  represented  by  R  and  0  ,  respectively.  The  radial  and 
azimuthal  magnetic  field  components  are  denoted  by  BR  and  BQm,  respectively.  The 
gradient  scale  length  at  the  magnetic  equator  is 


lb  *  IV<3'W3R!le  -o  ■  !  •  (3-1?> 

m 
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Consider  a  singly  ionized  monatomic  oxygen  ion  with  a  speed  of  1  x  10  km/s. 

The  appropriate  magnetic  field  strength  and  magnetic  field  gradient  scale  length 

4  -3  4 

for  R  =  3.5  x  10  km  at  the  magnetic  equator  are  1.87  x  10  G  and  1.17  x  10  km, 

respectively.  The  corresponding  maximum  possible  ion-gyroradius,  v/P,  is  8.9  x 
2 

10  km  and  so 


v/P  =  0.08  Lg  .  (3.13) 

Hence,  from  Eq.  (3.1)  it  follows  that  nonadiabaticity  can  preclude  the  long-term 
confinement  of  singly  ionized  monatomic  oxygen  ions  with  a  speed  of  1  x  10J  km/s  on 
the  magnetic  flux  surface  for  R  =  3.5  x  10^  km  at  the  magnetic  equator.  In  fact, 

p 

Eq.  (3.1)  implies  that  ion  speeds  smaller  than  5.2  x  10  km/s  are  probably  neces¬ 
sary  for  long-term  confinement  on  the  flux  surface  being  considered. 

Equations  (3.1),  (3.11b),  and  (3.12)  indicate  that  singly  ionized  ions  with 
an  atomic  mass  of  a  can  be  confined  on  flux  surfaces  with  a  radius  of  R-_(km)  at  the 
magnetic  equator  for  speeds  v(km/s)  satisfying  the  criterion 
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(3.14) 


,  .  1-0  x  1013 
v  <  2 
eq 

Equation  (3.14)  shows  that  the  ability  to  confine  energetic  ions  with  a  particular 

speed  is  decreased  by  higher  altitude  and  larger  atomic  mass.  For  example, 

confinement  of  singly  ionized  monatomic  oxygen  ions  with  speeds  of  1  x  103  and  2  x 

103  km/s  are  limited  to  magnetic  flux  surfaces  with  R  less  than  -2.5  x  104  km 
4  eq 

and  1.8  x  10  km,  respectively. 

From  the  standpoint  of  particle  kinetic  energy,  W(eV),  Eq.  (3.14)  implies 

that 


..  .  5.4  x  10 
W  <  — ■  4 


23 


(3.15) 


cxR 


eq 


Hence,  for  a  specified  kinetic  energy,  ions  with  a  lower  atomic  mass  can  be 
confined  to  a  higher  altitude  than  those  of  larger  atomic  mass.  For  example,  singly 
ionized  iron  ions  with  energy  below  -2.5  x  104  eV  can  be  confined  to  the  flux 

4 

surface  with  R  =  2.5  x  10  km,  while  the  corresponding  limit  for  singly  ionized 
eq  . 

mcnato'ric  oxygen  ions  is  -8.6  x  10  eV. 


3.5  SUMMARY  AND  CONCLUDING  REMARKS 

The  magnetic  moment  is  an  adiabatic  invariant  for  particle  motion  in  a 

3  1-34 

magnetic  field  geometry  with  sufficiently  weak  magnetic  field  gradients. 
Particles  with  sufficiently  large  speeds  can  have  trajectories  which  do  not 
preserve  the  constancy  of  the  magnetic  moment.  This  nonadiabatic  behavior  affects 
the  high-altitude  nuclear  environment  in  two  ways.  First,  the  magnetic  field 
geometry  can  be  significantly  altered  by  a  high-altitude  nuclear  burst  during 
early-times.  Nonadiabaticity  can  cause  pitch  angle  scattering  of  ions  at  points  of 
maximum  (minimum)  magnetic  field  strength  in  this  modified  geometry.  This  process 
supports  the  premise  of  minimum  outer  scale  length  in  the  high-altitude  nuclear 
environment  being  at  least  as  large  as  the  energetic  ion-gyroradii.  Second,  ions 
injected  into  the  magnetosphere  following  a  high-altitude  nuclear  burst  can  be 
confined  In  the  magnetosphere  only  If  their  kinetic  energy  is  below  a  bound  deter¬ 
mined  by  the  onset  of  nonadiabaticity.  This  result  implies  that  the  characteriza¬ 
tion  of  the  ionic  population  in  the  magnetosphere  and  the  effect  of  the  ions  on  the 
ionosphere  is  influenced  by  the  nonadiabatic  scattering  associated  with  the 
magnetospheric  magnetic  field. 
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